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ABSTRACT 

We extract a sample of disc stars within 200 pc of the Sun from the RAVE and SDSS 
surveys. Distances are estimated photometricahy and proper motions are from ground- 
based data. We show that the velocity-space substructure first revealed in the Geneva- 
Copenhagen sample is also present in this completely independent sample. We also 
evaluate action-angle variables for these stars and show that the Hyades stream stars 
in these data are again characteristic of having been scattered at a Lindblad resonance. 
Unfortunately, analysis of such local samples can determine neither whether it is an 
inner or an outer Lindblad resonance, nor the multiplicity of the pattern. 

Key words: Galaxy: disc - Galaxy: kinematics and dynamics - solar neighbourhood 
- Galaxy: structure - galaxies: spiral - stars: distances 



THE LOCAL DISTRIBUTION OF STELLAR 
VELOCITIES 



Analysis of the Hipparcos data by iDehnenI (|l998h re- 
vealed that the stellar velocity distribution in the so- 
lar neighbourhood manifested significant substructure. 
The heroic Geneva-Copenhagen Survey (he reafter GCS 
iNordstrom et al\ |2004| : iHolmberg et al\ I2OO9I ) followed up 
with radial velocity measurements of 14 139 nearby F and 
G dwarf stars, which confirmed and strengthened Dehnen's 
conclusion. 

The Geneva-Copenhagen survey was constructed so as 
to avoid most of the selection biases that went into the full 
Hipparcos sample. Aside from a concentration of 112 stars in 
the Hyades cluster, the distribution of sample stars over the 
sky is remarkably uniform, with a slightly higher density in 
the declination range south of 5 = —26°. The large majority 
of stars are within 200 pc of the Sun and the sample within 
40 pc is believed to be nearly complete. 

The Cartesian heliocentric velocity components of stars 
near the Sun in Galactic coordinates are U, V and W, with 
U being directed towards the Galactic centre, V being in 
the direction of Galactic rotation, and W towards the north 
Galactic pole. The velocity substructure in the GCS is par- 
ticularly evident in the U—V plane, where the distribution is 
broken into a number of substantial streams, with no under- 
lying smooth component. Numerous studies ( Famaev et al. 



l201lD have shown that the streams are both too substantial 
and chemically in homogeneous to be dissolved star clusters 
(e.Q. lEggen| [l99^. 

Dynamical evolution is the most likely source of 
the features, which have been modelled extensively. 
iDe Simone et al\ (|2004l ) suggest that the entire velocity dis- 
tribution arises fr om a succes s ion o f transient spiral per- 
turbations, while iHelmi et al\ (|2006D attribute some sub- 
structure to minor accretion events. Individual features 
have been modelled as responses to the bar and /or vari- 
ous assumed spiral perturbations within the disc (iDehnen 
2000| : lQuillenll2003l. jQuillen fc Minchevll2005l: IChakrabartvl 



lying smooth component. Numerous studies (|l 
l2007l : iBensbv et af.l2007l : IBovv fc Hog"3l2010l : E 



'ompeia et al. 
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20071 : lAntoia et aZ.ll2009l: Minchev et al\\20m . The anal- 
ysis of Quillen fc Minchev l|2005h . which was extended by 
[Pompcia et al. (201^ identified the inner ultra-harmonic 
resonance of an assumed 2-arm spiral as a possible cause of 
both the Hyades and Sirius streams. 

By contrast, ISellwoodI (|2010l . hereafter Paper I) pre- 
sented an analysis that did not need to assume a form for 
the perturbation. Using action-angle coordinates, he showed 
that the stars of the Hyades stream were both concentrated 
along a resonance line in action space and grouped in a com- 
bination of angle coordinates (possibly) ind icative of a recent 
inner Lindblad resonance (hereafter ILR). iMcMiliaiil (1201 ll ) 
has confirmed most of his analysis, and concurs that a Lind- 
blad resonance is responsible but, unfortunately, was able to 
show that subtle selection effects in such local data imply 
that the distribution could also be consistent with trapping 
in an outer Lindblad resonance (OLR). We discuss selection 
effects and McMillan's analysis in the Appendix. 

While the GCS sample offered compelling evidence for 
a Lindblad resonance, it is desirable to attempt to confirm 
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the result from other indep endent surveys. The R adial Ve- 
locity Experiment fRAVE. Isteinmetz et aZ.ll2006l ). a large 
spectroscopic survey of stars in the southern sky, plans to 
measure the heliocentric radial velocities and stellar param- 
eters for about a million stars in the apparent magnitude 
range 9 < mi < 12; the fi rst 21 121 were ma de available 
in the second data release (jZwitter et aZ.llioOsh . The typi- 
cal uncertainty in the radial velocity is < 2 km s~^, but 
the distance to most stars has to be judged photometrically 
and most proper motions are from ground-based data. Thus 
the three phase space coordinates for each star are of much 
lower quality than are those in the Geneva-Copenhagen sur- 
vey, although this weakness will, when the survey is com- 
plete, be compensated by a much larger sampl e size. The 
huge northern Slo an Digital Sky Sur vey (SDSS, lYork et all 
|200C| ) and Segue2 (lYannv et a/.li2009l l surveys are complete, 
but sample fainter stars (the magnitude range for Segue2 
was 14.0 < g < 20.3) that are therefore generally more 
distant than are the RAVE stars. The recent ly-released M- 
dwarf sample of SDSS stars (|West et ttZ.ll20lil ) substantially 
increases the number of stars with estimated distances and 
kinematics within the neighbourhood of the Sun. 



2 SAMPLE SELECTION 
2.1 RAVE stars 

We have downloaded the on-line table of the second data re- 
lease from the RAVE website and selected a subset of stars 
for analysisQ We estimate dis tances to these stars b y fitting 
to the Yonsei-Yale isochrones (jPemaraue et al using a 

method related to that described bv lBreddels et aU l|2010h . 
We adopt many of their selection criteria: we require the 
spectral signal-to-noise parameter > 20 with a blank spec- 
tral warning flag field; the parameters [M/H], log(gr), and 
Teff to be determined; and the stars to have J & Ka m agni- 
tudes from two-micron aU sky survey ('2MASS[s| [200^ with 
no warning flags about the identiflcation of the star or the 
2MASS photometry. Unlike those authors, however, we have 
kept stars with b < 25° on the grounds that extinction for 
the nearby stars that interest us will not be large enough 
in the near IR to severely bias our distance estimates. As 
we wish to select nearby main-sequence stars that are mem- 
bers of the disc population, we also eliminate stars with 
log(g) < 4, Toff > 10* K, and with \vr\ > 80 km s"\ 

We estimate the absolute J magnitude of each selected 
star by matching the estimated [Fe/H], log(g), Ten, and J-Ks 
colour to values in the isochrone tables for stars of all ages 
and all values of [a/Fe], rejecting a few more stars for which 
the best match > 6. We consider the closest match in 
the tables to the given input parameters to yield the best 
estimate of the absolute magnitude from which we estimate 
a photometric distance using the apparent J-band magni- 
tude0 Our resulting sample contains 7 384 stars. We save 
the values of [Fe/H], log(,g), T^b, and J-Ks colour of the 



^ The third release oc curred while this paper was under review. 

IZwitter et al. I l l20ld) describe a similar method, but define a 
"most likely" estimate of the absolute magnitude that differs from 
our "best" estimate. The difference is likely to be well within the 
uncertainties for the main-sequence stars considered here. 
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Figure 1. Histograms of distances and (7-velocity uncertainties 
for the selected sample of 5145 RAVE stars. 



closest matching model star in the isochrone table; Monte- 
Carlo variation of the stellar parameters about this saved set 
of values suggests that distances have a relative precision of 
30% - 50%, with some larger uncertainties. 

We use the proper motions in equat orial c oordinates 
tabulated in RAVE, mostly from Tycho-2 (|H0g et al..200d ). 
which we then combine with the radial velocity and position 
to determine the heliocentric velocity in Galactic compo- 
nents U, V SzW l|johnson fc Soderblomlll987l : iPiatek et all 
l2002h . We estimate uncertainties in these velocities from 
5000 Monte Carlo re-selections of all the stellar parameters 
that affect the distance estimate, adopting cr(J) — 0.03 mag, 
a(J-K,) = 0.042 mag, cr(rcff) = 300 K, a(lo£q) = 
0.3 dex, and a([Fe/H]) = 0.25 dex ijBreddels et aimOld ). 
as well as the tabulated radial velocity and proper motion 
uncertainties. 

In order to select nearby thin disc stars, we further re- 
strict the sample to stars whose best estimate of the dis- 
tance is within 500 pc and retain only those having an 
energy of vertical motion about the Galactic mid-plane, 
= 0.5{z'^u^ + W^) < 392 (km s~^)^, with the vertical fre - 
quency v = 0.07 km s"^ pc"^ l|Binnev fc Tremaind lioosl ) . 
giving them a maximum vertical excursion of ±400 pc. This 
latter restriction rejects most thick disk and halo stars, as 
well as a few thin disc stars. 

Fig. H] shows histograms of distances and of velocity 
uncertainties for the 5, 145 remaining stars. The mode of the 
distance distribution is 200 pc from the Sun. The estimated 
uncertainties in U are generally < 10 km s"'^, with a tail up 
to values four times larger. 



2.2 SDSS and Segue2 stars 

We also select ed stars from the 
(jAbazaiian et al\ hoOQ) with 4 < 
3000 < Toff < 10* K, and with 



DR7 of SDSS 
log(g) < 5, 
< 60 km s~^, as 



estimated by the Sloan spectral parameters pipeline 
(|Lee et ttZ.ll2008l , SSPP), which yielded over 100 000 candi- 
date stars. As for RAVE, we wish to estimate photometric 

dist ances to these star s. 

lJuric et al\ l|2008h and llvezic et al\ (|2008l ). respectively, 
propose formulae for estimating the absolute magnitudes, 
Mr, from the r — i colour, or from the g — i colour and [Fe/H] 
but, as for the RAVE stars, we have preferred to estimate 
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Figure 2. Blue points show the estimated Mr magnitude as 
a function of colo ur and metallicity using the formulae from 
llvezic et al. Red points show the same quantity obtained 

from isochrone fitting, while green triangles show values for which 
no acceptable isochrone fits the star. See text for discussion. 

distances by fitting tlie full stellar para meters to isochrone 
tables. The Dartmouth isochrone tables jPotter et aZ.|[200^ ) 
give absolute magnitudes in Sloan colour bands. We inter- 
polate for [Fe/H], for all ages and with [a/Fe]=0 only, to 
estimate Mi by fitting r — i colours, together with SSPP 
estimates of [Fe /H], log((?), and Tcfi. As recommended by 
I Lee et al\ l|2008l ). we increased the uncertainties in the es- 
timated stellar parameters to a(T'eff) = 157 K, cr(log g) = 
0.28 dex, and cr([Fe/H]) = 0.24 dex. We find 38 690 stars 
whose best estimate of distance modulus in the i-band places 
them within 2 kpc. 

Fig.[2]illustrates the superiority of isochrone fitting over 
relying on a single colour index. The blue points show the 
estimated Mr from the fifth-order polynomial function of 
the g — i colour, broadened by a spread g iven by a quadratic 
expression in [Fe/H], as recommended bv llvezid et al\ (|2008l . 
eqs. A2 & A7). The red points show the isochrone-fitted Mr 
magnitude obtained as described above; note that the SSPP 
does not return parameters for stars having Tefi < 4 500, 
which is the reason for the absence of credible main sequence 
stars redward of g — i ~ 1.7. For the large majority of stars, 
the blue and red points fill the same region of the Figure, 
and are therefore consistent; the two magnitude estimates 
differ by < 0.2 mag for 89% of the stars. But, as is physi- 
cally reasonable, isochrone fitting broadens the distribution 
near the main sequence turn off, although the spread may 
be spuriously enhanced by errors in the pipeline parameters. 




Figure 3. Histograms of distances and (7-velocity uncertainties 
for the selected sample of 629 SDSS stars. 
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Figure 4. Histograms of distances and (7-velocity uncertainties 
for the selected sample of 10 669 M-dwarf stars. 



especially in \og{g). The principal advantage of isochrone fit- 
ting, however, is that it flags stars with inconsistent param- 
eters. The green triangles show 938 stars that had > 6, 
based on the above uncertainties, for the smallest difference 
between the estimated r — i colour, Tes, and log(gi) of any 
model star in the Dartmouth isochrone tables. A large 
generally arises for a star that has a colour (we used r — i) 
estimate inconsistent with the pipeline-estimated Tcff, al- 
though the estimated log (3) may be inconsistent in other 
cases. We therefore adopt the absolute magnitudes given by 
the isochrone fits to estimate distances, discarding the stars 
having a large . 

We adopted the distance modulus in i, and limited the 
distance and z motions as for RAVE stars. These restric- 
tions reduced the sample from the SSPP to just 629 nearby 
disc stars. This disappointing immber results from the ap- 
parent magnitude limit for SDSS spectroscopy of rrii > 14 
coupled with the elimination of all M-dwarfs from the sam- 
ple through the low-temperature restriction of the SSPP; 
the distance histogram shown in Fig. [3] indicates that no 
star from the SSPP is closer than 200 pc. Because uncer- 
tainties in the stellar parameters are somewhat smaller than 
for RAVE stars, the estimated uncertainties in distance are 
slightly smaller, which translates into smaller velocity un- 
certainties. 
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Figure 5. Top right: the sky distribution in Galactic coordinates of the selected RAVE, blue, and SDSS M-dwarf, red, stars within 
200 pc of the Sun. Top left: the Galactic velocity components of the selected stars with the same colour coding as in the top panel. 
Bottom left: contours of the density of the combined sample in velocity space. The shapes of the contours are quite similar to those for 
the Geneva-Copenhagen sample reproduced from Paper I in the bottom right panel, where the letters indicate the approximate locations 
of the principal "star streams": A — Hyades, B - Sirius, C — Hercules, and D — Pleiades. 



2.3 M-dwarf sample 

IWest et al] |2oi3) provided a catalogue of 70 841 M-dwarfs 
from SDSS DR7 that were omitted from the SSPP. Their 
table provides a photometric estimate of the distance to each 
star, as well as a radial v elocity and proper motio n from 
USNO-B/SDSS catalogue (|Munn et aLll2004 |2008| ). They 
suggest that distance uncertainties are typically about 20% 
and uncertainties in radial velocity are 7-10 km s~^. 

As IWest et al\ (|201ll ) recommend, we selected stars 
with the "goodPM" and "goodPhot" flags set to 'true', and 
the "WDM" flag set to 'false' to eliminate possible binaries 
with a white dwarf companion, which reduces the sample 
to 39 151 stars. We further excluded stars having no radial 
velocity as well as those with no distance estimate or for 
which the estimated distance exceeded 500 pc. Finally, we 
also eliminated stars for which any component of the he- 
liocentric velocity exceeded 80 km s"'^ and those having a 



vertical energy that would take them farther than 400 pc 
from the disc mid-plane, leaving us with a final sample of 
10,669 stars. 

We estimated uncertainties in the velocity components, 
U , V &Z. W hj combining the 20% distance uncertainty, a 
10 km s~^ uncertainty in the radial velocity, and the proper 
motion uncertainties. The distribution of distances and U 
velocity errors is shown in Fig. |4l many of these intrinsically 
faint stars lie within 200 pc and, again, velocity uncertainties 
are typically 10 km s~^. 



2.4 Combined sample 

The three surveys yielded a total sample of some 16 443 stars 
within 500 pc. Since few stars in these surveys are close to 
zero Galactic latitude, the in-plane velocity components U 
and V are mostly determined by the proper motion, with 



© 2010 RAS, MNRAS 000, 



Velocity- space substructure 5 



the more-precise radial velocity making a smaller contribu- 
tion. For this reason, we confine all the following analysis to 
the 6, 769 stars closer than 200 pc to the Sun, which have 
correspondingly smaller uncertainties in these two velocity 
components. We refer to this as our final "combined sam- 
ple". 

Note that the combined sample contains just four stars 
from the RAVE survey, and no stars from the SDSS, that 
were also in the GCS. The reason for this tiny overlap is that 
the GCS was limited to F and G dwarf stars brighter than 
JTivis = 8.6, and the only stars from SDSS in the combined 
sample are the M-dwarfs. Thus the two samples are truly 
independent; even the radial velocities of the four stars in 
common were remeasured by RAVE. 

The top right panel of Fig. [S] shows that the combined 
sample covers most of the sky, except for the Galactic plane, 
although far from uniformly. The top left panel shows the 
best fit velocity components in Galactic coordinates. Al- 
though velocity uncertainties are several times larger than 
for the GCS, the general appearance of the distribution in 
velocity space is similar. 

The bottom left panel of Fig. [S] contours the density in 
the space of these two velocity components, using linearly- 
spaced contour levels. The bottom right panel reproduces 
the plot, with logarithmically spaced contours, constructed 
from the GCS sample that was shown as Fig. 2 in Paper 
I. Many of the features in the ([/, V)-p\a,ne that stand out 
in the GCS have counterparts in this completely indepen- 
dent sample, although their velocity locations do not match 
perfectly. The strongest feature is the Hyades stream, but 
the bottom left panel has clear hints of the Pleiades, Sirius, 
and even Hercules streams. Note that errors in our distance 
estimates will give rise to correlated errors in the U and 
V velocity components, since they are derived mostly from 
proper motions. However, the misplacements of points in the 
{U, 1/)-plane are in directions that differ for stars in differ- 
ent parts of the sky, leading only to a general blurring of 
feat ures. 

iKlement et~al. I l|2008l . 2011) present a more sophisti- 
cated analysis for possible streams am ong nearby stars in 
the fir st two RAVE data releases, and iKaratas fc KlementI 
l|201ll ) have recently confirmed many of the features in the 
{U, V) plane from an analysis of RAVE stars selected to have 
smaller velocity uncertainties. 



3 ANALYSIS 

We compute action-angle variables (|Binnev fc Tremaind 
|2008| ) for the in-plane motions of the selected stars, as de- 
scribed in Paper I, using our best estimates of the full phase- 
space coordinates for each star, co rrected for the motio n 
of the local standard of rest (LSR ISchonrich et al\ |2010| ) . 
Monte Carlo simulation using the uncertainties in the in- 
put data indicated that the median action uncertainties are 
(j{Jr) = 0.0032 and o{J^) — 0.020, while the median angle 
uncertainties are a{wii) = 0.18 and <j{w^) — 0.024 radians. 
The dimensionless actions are scaled by Lzfi = RoVo, where 
Ro is the solar distance from the Galactic centre and Vb is 
the orbital speed of the LSR, in order to make them inde- 
pendent of these two uncertain quantities. The upper panel 
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Figure 6. The upper panel shows the distribution of actions esti- 
mated for the 6 769 stars of the combined sample. The parabolic 
lower boundary is a selection effect, since only stars with larger 
radial oscillations can visit the solar neighbourhood if their guid- 
ing centre radii differ from Rq. 

The lower panel shows the mean density of stars in the up- 
per panel along resonance lines for the ILR of m = 2 distur- 
bances having a range of pattern speeds. The shaded region shows 
the 99% confidence range from randomly resampled coordinates, 
while the solid line shows the same quantity from the selected 
sample. The equivalent plot to the lower panel shows the same 
significant feature, except at a quite different frequency, when we 
assume trapping at either an ILR or an OLR of a 2-, 3- or 4- fold 
symmetric pattern. 

of Fig. [6] shows the distribution of actions for the stars in 
the combined sample. 

Stars in resonance with a weak perturbation having m- 
fold rotational symmetry and which rotates at the angular 
rate ilp have unforced frequencies that obey the relation 
mSl^ -I- irin = mQ,p, where I = at corotation, and / = ^1 
at the ILR and OLR respectively. This frequency condition 
defines a line in action space that slopes to smaller for 
increasing Jr for all {l,m). 

The lower panel shows the density of stars along ILR 
(Z = — 1) lines for m = 2 disturbances having a range of pat- 
tern speeds. As in Paper I, we estimate the significance by 
comparison with equal size samples of pseudo-stars but with 
radial coordinates and in-plane velocities chosen indepen- 
dently in a bootstrap fashion from the distributions of these 
three variables. The shaded area encloses 99% of the results 
from these randomly resampled coordinates, while the solid 
line shows the same quantity from the selected sample. The 
peak value of the solid curve is thus statistically highly sig- 
nificant and occurs for f2p = 0.294 (in units of Vo/i?o), in 
excellent agreement with the value (0.296) obtained from 
the completely independent GCS sample. 

The excess at this frequency is caused by the overdense 
feature visible in the upper panel slanting upward with neg- 



© 2010 RAS, MNRAS 000, [TH9] 



6 Hahn, Sellwood and Pryor 



0,4 - 

>^ 

£ 0.3 - 








Figure 7. Distributions of radial pliase angles ui^ (upper panel) 
and of orbital phase angles (lower panel) for the combined 
sample. 



ative slope from near the point (1,0). Unfortunately, this fea- 
ture cannot be uniquely attributed to an ILR of an m = 2 
pattern, since the loci of corotation and of both Lindblad 
resonances for m J5 2 also have similar steep negative slopes 
in action space, as shown in the lower panel of fig. 5 of Paper 
I. Thus tests for trapping at any resonance appear very sim- 
ilar to that shown in the lower panel of Fig. 6, although the 
almost equally significant excess when other resonances are 
considered clearly must imply quite different pattern speeds. 

Stars that are scattered at a resonance, on the other 
hand, change both actions in such a way that the Jacobi con- 
stant is conserved, which generally shifts a star in a direction 
in action space that is not parallel to the resonance locus. 
Only at the ILR does the scattering line in action-space have 
similar, though not identical, slope to the resonance locus; 
scattering lines at corotation are horizontal while they have 
positive slope at an OLR. Tests for an excess of stars along 
lines of positive slope in action space, that would correspond 
to scattering at an OLR, confirmed that there are no signif- 
icant features with this slope, as was also true for the GCS 
sample in Paper I. 

Figure[7]shows the distributions of the two angles conju- 
gate to the actions. (The physical meaning of these variables 
is explained in the Appendix.) The orbital phase distribu- 
tion, w^, is narrow because the guiding centres of all stars 
in the sample are close to the Sun's azimuth in the Galaxy, 
which was arbitrarily chosen to lie at = 0. The distribu- 
tion of radial phases, wr, is non-uniform also, reflecting the 
substructure in phase space (Fig. [5]). 

As explained more fully in Paper I, a new peak that 
appears in any of the distributions of simple linear combi- 
nations of these phase angles may indicate a group of stars 
trapped in, or recently scattered by, a resonance. The ap- 
propriate combination is mwtf, -\- Iwr . This test is insensitive 
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Figure 8. Tests of the distributions of angles for trapping at 
an outer Lindblad resonance of an m = 2 disturbance (upper 
panel) and respectively of an ILR for m = 2, m = 3 and m = 4 
disturbances (lower three panels). 



to corotation, where I = 0, but a new concentration of stars 
at some value of one of these combinations with \l\ = 1 is 
an indicator of a Lindblad resonance. McMillan! (|201ll ) has 
shown that selection effects in any local sample, together 
with a small amount of scatter about the expected constant 
value of mw^ + Iwr, conspire to frustrate this test. In the 
light of his finding, our more limited objective here is to 
show that the features that appeared in the tests in Paper 
I have their counterparts in the present sample. 

Fig. [8] shows the distributions of four combinations of 
the angle variables of stars in the combined sample. The 
top panel shows the case for an OLR for m = 2, while the 
lower three panels are for various ILRs. The overall shapes 
of the distributions differ in detail from those found for GCS 
stars in Paper I. However, Sellwood identified clear peaks 
in the ILR cases for m = 2 and m = 3 (his Fig. 7), and 
for m = 4 in Fig. IA3I of this paper. These features, at the 
abscissae marked by arrows, can also be identified in Fig. [8] 
in this independent sample. Although the significance of the 
peaks at the indicated abscissae is not high, they do lie at 
the same phases as those in the GCS sample. The peak near 
zero, that becomes more prominent as m is increased, is a 
selection effect that is explained in the Appendix. A peak 
is also visible in the OLR distribution, in the top panel, 
near 2w0 + wr ^ —1.8 that also has a counterpart at the 
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indicated position in the corresponding distribution from the 
GCS sample in Paper I, although in that case the difference 
between the distributions of 2w,j, + wr and of wr alone was 
less pronounced. 

The peaks in these new data are far from compelling, 
but they provide a valuable confirmation of the far more 
significant peaks that were found from the higher quality 
GCS sample. Measurement errors must always smooth away 
features on the scale of the uncertainty and the above esti- 
mates of the uncertainties suggest broadening on the scale 
of (J ~ 0.2 radians. In particular, it is reassuring that the 
most significant peak in Fig. [5] lies at the same frequency, 
and the peaks in the lower three panels of Fig. [S] lie at the 
same phases, as in the corresponding Figures from the GCS 
sample. Taken together, this is strong evidence in support 
of a Lindblad resonance from this completely independent 
sample. 



dense \ow-Ju region up to higher Ju where the overdensity 
stands out. Note that if the cause is a spiral (other types 
of disturbance could be responsible), the ILR of a bisym- 
metric spiral seems unlikely, since it would place corotation 
unreasonably far out in the disk, but an m = 4, or perhaps 
even an m = 3, spiral would seem more likely, as noted in 
Paper I. This speculat ion could be tested by data from Gaia 
(|Perrvman et a t 2001), which will obtain phase space infor- 
mation over a much more extensive region of the Galaxy. 
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4 CONCLUSIONS 



REFERENCES 



We have show n that the velo city-space sub structure revealed 
by H i pvarcos llDehnenI 199^) and the GCS (|Nordstrom et al\ 
12004 iHolmberg et all |2009| '1 is also present among the 
nearby (d < 200 pc) stars in the RAVE and SDSS/Segue2 
samples. We find that the velocity space substructure closely 
resembles that found for the independent sample of GCS 
stars even though the present sample is barely half the size 
and velocity uncertainties are several times larger; most of 
the significant star "streams" of the GCS have counterparts 
in the new sample, although the precise velocities of the 
streams do not match perfectly. 

Analysis of the action variables constructed from our 
best estimates of the full phase space coordinates of each 
star in the sample reveals an excess of stars along a reso- 
nance scattering trajectory in actio n-space at a sim ilar fre- 
quency to that found in Paper I ISellwoodl [2OI0I ) for the 
GCS sample. While still statistically highly significant, the 
feature in Fig. 6 stands out less clearly than in Paper I be- 
cause the sample is smaller and uncertainties are larger. The 
evidence for resonant trapping in Fig [8] is also weaker than 
that found for the GCS sample, but again reveals peaks at 
the same phases. We therefore consider this sample of stars 
to support the evidence for a Lindblad resonance found in 
Paper I, although here again the data do not constrain the 
multiplicity of the pattern. 

While it is unfortunate that these data do not rule out 
trapping at an O LR as the cause of the Hyades stream, it is 
worth noting that lDehnenI (|2000t ) examined the effect on the 
local phase space due to the OLR of the bar and did not find 
any Hyades stream-like features. Since the Hyades stream 
stars form the tongue in action-space that stands out Fig. 3 
of Paper I and our Fig. (6] it seems to us far more like l y that 
the stream was created by an ILR. ISellwoodI l|l994 hood ) 
reported that features that extend upward towards smaller 
Jtf, (= Lz) for increasing Jr (or -Erand) are created by ILRs 
in his simulations. Since scattering at an OLR moves stars in 
action space in a direction that is roughly perpendicular to 
the resonance locus, as discussed in §3 above, stars do not 
move far before leaving the resonance. Only in the case of an 
ILR do stars stay close to resonance as they get pushed by 
the disturbance, thereby allowing stars to be moved from the 
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Figure Al. Four possible orbits that pass through the vicinity of 
the Sun (asterisk), drawn in the frame that rotates at the angular 
rate of the LSR. The sense of rotation of the Galaxy is clockwise, 
in line with the usual convention. The sample volume is the lightly 
shaded region, and the dotted line marks the solar circle. 
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APPENDIX A: SELECTION EFFECTS 

Both the present sample and the GCS sample of Paper I 
were confined to stars that happen to lie near the Sun at 
present. Since stars move on eccentric orbits, most stars in 
the sample are merely visiting the solar vicinity from farther 
afield. 

Fig. I All shows four possible orbits of stars near the Sun 
that all happen to pass through the (shaded) survey vol- 
ume. The orbits are drawn in the frame of the LSR and have 
sufficiently small radial excursions to approximate Lindblad 



Figure A2. The distribution of angles for GCS stars. Stars with 
radial action > 0.05Lz^o a-re excluded. Note the difference in scales 
between the two axes. 



epicycles. All four orbits have the same Ju (epicycle size); 
those labelled 2 & 4, with guiding centres (plus symbols) 
lying on the solar circle, have the same {— Lz) as does 
the Sun, while orbits 1 & 3 have, respectively, greater and 
smaller values of J^. Orbits 2 & 4 are closed ellipses, be- 
cause their guiding centres move at the same angular rate 
as the frame, but the guiding centres of the other two el- 
lipses, which are marked at the point of the star's closest 
approach to the Sun, drift relative to the LSR. For clarity, 
the orbits shown do not loop around the Sun so that the 
vectors, which show peculiar velocities relative to the LSR, 
appear to correlate with the direction of the star from the 
solar position. In reality, stars anywhere in the survey region 
can have peculiar velocities in any direction. 

The value of ui^ is the instantaneous G alactic azimuth 
of the guiding centre, which ISellwoodI l|2010l . Paper I) chose 
to be zero on the line connecting the Sun to the Galactic 
centre. Thus > for stars ahead of the Sun. The wr 
variable is the phase of the star around its epicycle, which 
increases with t ime in the sense shown by the arrows. Again 
ISellwoodI HoiO) adopted wr = for a star at its apocentre 
and consequently wr = ±7r for a star at its pericentre. Note 
that action-angle variables, unlike Lindblad's epicycles, can 
also be used for more eccentric orbits, for which the guiding 
centres generally orbit more slowly than the circular speed 
and the radial oscillation is anharmonic. 

Fig. 4 of Paper I showed the separate distributions of wr 
and for stars in the GCS sample having Jr < 0.05L2,o- 
Here, Fig. IA2l shows the joint distribution of the same stars 
in the space of both angles. The azimuthal phase distribution 
is confined to |w^| < 0.4 reflecting the limited spread in 
Galactic azimuths of the guiding centres for stars that pass 
through the survey volume. 

The distribution of values for stars having wr = 0, 
such as orbit 3 in Fig. lAll or wr = zbvr (orbit 1) is very 
narrow because their guiding centres lie along the radius 
vector from the Galactic centre to the Sun for all values of 
Jr. Stars in the first quadrant of Fig. IA2l are inward moving 
stars (because wr > 0) with guiding centres ahead of the 
Sun {w^ > 0) , such as orbit 4. Conversely, stars in the third 
quadrant (ui^ < 0) are outward moving stars having wr < 
{e.g. orbit 2). Thus the selection criteria preclude stars from 
lying in the second and fourth quadrants. The spread of w^ 
values for wr ~ ±7r/2 refiects the spread in the sizes of 
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Figure A3. The distribution of 4i«^ — wn for GCS stars, whicli 
was not shown in Paper I. As there, stars with radial action > 
O.OSL^ a-re excluded. 



the epicycles of the stars in the sample, from which more 
ecce ntric orbi t s were eliminated. 

ISellwoodI l|2010l ) argued that a concentration of stars 
having a nearly constant value of mw^ + Ivor would be an 
indicator of trapping in a resonance for the selected values of 
I and m. Thus he searched for an excess of stars lying along 
lines of fixed slope with all possible intercepts in this plot, 
and reported the results in his Fig. 7 of Paper I for the GCS 
sample. The lines for inner Lindblad resonances {I = — 1) 
have a positive slope in Fig. IA2I Their slopes decrease as m 
increases, causing them to become more closely parallel to 
the distribution near (we., w^) = (0, 0), giving rise to a peak 
near mw^ — wr — as m increases. The distributions shown 
in Fig. 7 of paper I for m = 2 and m = 3, and in Fig. IA3l of 
this paper, show the increasing prominence of this feature, 
which is purely an artefact of the sample selection. A similar 
feature appears near mw^ +wr = ivr when I — 1, since the 
line s have negative s lope in these cases. 

ISellwoodI (|2010l ) drew attention to the peaks that lay 
away from mw^ — wr = for the cases of inner Lindblad 
resonances. These peaks arise in part from the concentration 
of GCS stars near wr ~ —1.5, ~ —0.15 in Fig | A 2l but the 
clump is visibly extended, which iMcMillanI (|2011 ) described 
as having a triangular shape. Since the excess lies in the 
third quadrant, it gives rise to peaks in the distributions 
of mwtj, — Wr that shift closer to zero as m is increased. 
Obviously the same excess gives rise to peaks near mw^ + 
Wr > —2 for lines of negative slope in the cases of an OLR 
but, at least for the GCS sample, they are not as striking as 
those for the ILR cases. 

The case for an ILR made bv ISellwoodI (|2010l 'l did not 
rest solely on this slight difference, however. He attached 
far greater significance to the fact that stars lying along a 
line of positive slope extending through the clump were ex- 
actly those that lay along the resonant scattering tongue 
in action space, as shown in Fig. 8 of Paper I. As this was 
not the case for stars lying on lines of negative slope, Sell- 
wood concluded that an ILR was responsible. Unfortunately, 
iMcMillanI l|201l[ ) found that his conclusion was compromised 
by selection effects: if stars are not only in the vicinity of the 
Sun but are also trapped in a resonance, McMillan was able 
to show that both ILR and OLR models give rise to fea- 
tures in the angle distribution {i.e., Fig. lA2|) having positive 
slope. Thus, although scattering at a Lindblad resonance 
was clearly established, the test in Paper I did not enable 
the character of the resonance to be determined. 
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